Baryon Stopping and Charged Particle Distributions in 
Central Pb+Pb Collisions at 158 GeV per Nucleon 
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Net proton and negative hadron spectra for central Pb+Pb 
collisions at 158 GeV per nucleon at the CERN SPS were 
measured and compared to spectra from lighter systems. Net 
baryon distributions were derived from those of net protons, 
utilizing model calculations of isospin contributions as well 
as data and model calculations of strange baryon distribu- 
tions. Stopping (rapidity shift with respect to the beam) and 
mean transverse momentum (pr) of net baryons increase with 
system size. The rapidity density of negative hadrons scales 
with the number of participant nucleons for nuclear collisions, 
whereas their (pr) is independent of system size. The (pr) 
dependence upon particle mass and system size is consistent 
with larger transverse flow velocity at midrapidity for Pb+Pb 
compared to S+S central collisions. 

PACS numbers: 25.75-q,25.75.Dw 



Lattice QCD predicts that strongly interacting mat- 
ter at an energy density greater than 1-2 GeV/fm 3 at- 
tains a deconfined and approximately chirally restored 
state known as the quark-gluon plasma (for an overview, 
see Q). This state of matter existed in the early uni- 
verse, and it may influence the dynamics of rotating neu- 
tron stars The collision of nuclei at ultrarelativistic 
energies offers the possibility in the laboratory of creat- 
ing strongly interacting matter at sufficiently high energy 
density to form a quark gluon plasma || . Hadronic spec- 
tra from these reactions reflect the dynamics of the hot 
and dense zone formed in the collision. The baryon den- 
sity, established early in the reaction, is an important 
factor governing the evolution of the system Q. Com- 
parison of model predictions with measured rapidity and 
transverse momentum distributions and correlation func- 
tions constrains the possible dynamical scenarios of the 
reaction |5| , such as those for longitudinal and transverse 
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flow p] . In addition, the mechanism by which the incom- 
ing nucleons lose momentum during the collision (baryon 
stopping 0) is an important theoretical problem [p|^0|, 
and the measurement of baryon stopping in heavy ion 
collisions provides essential data on this question. 

In this letter, we present measurements by the NA49 
collaboration of rapidity and transverse momentum dis- 
tributions of participating baryons and negative hadrons 
over a large fraction of phase space for central Pb+Pb 
collisions at 158 GeV per nucleon (yijvw = 17-2 GcV). 
Hadronic spectra from S+S collisions at 200 GeV per 
nucleon pl|-p^[ and N+N (nucleon-nucleon) collisions at 
200 and 400 GcV Jl4|,[l5| serve as important references, 
helping to identify effects that depart from those ex- 
pected from the linear superposition of many N+N col- 
lisions. 

The NA49 apparatus is described in A beam of 
208 Pb ions at 158 GeV per nucleon (yi a b = 5.8) from the 
CERN SPS struck a 224 mg/cm 2 thick natural Pb tar- 
get. The 5% most central collisions were selected by mea- 
surement of the forward-going energy in the phase space 
occupied by the projectile spectator nucleons. The reac- 
tion products passed through a dipole magnetic field, and 
tracks from charged particles used for this analysis were 
measured in two large Main Time Projection Chambers 
(MTPCs) placed downstream of the magnets on either 
side of the beam axis. Over 1000 tracks were measured 
in each event. Identification of protons utilized the spe- 
cific ionization ((dE/dx)) of the gas of the MTPC @. 
A relative (dE/dx) resolution of 6% was achieved. 

The evolution of the incoming baryons (baryon stop- 
ping) was studied through measurement of the differ- 
ence of the proton and antiproton distributions (net pro- 
tons, denoted p—p) to eliminate the effect of baryon- 
antibaryon pair production. Meson production was stud- 
ied through the yield of negative hadrons (/i~), compris- 
ing primarily ir~ , with an admixture of K~ and p. 5xl0 4 
central events were used for the p—p analysis and 7xl0 3 
central events were used for the h~ analysis. 

The net proton yield was determined by an identifica- 
tion technique designed to minimize the systematic er- 
rors over wide acceptance Q. For each phase space 
bin (y,pr)i the distribution of (dE/dx) from negatively 
charged particles was subtracted from that of positively 
charged particles, resulting in a distribution that con- 
tained a peak with positive amplitude (more positive 
than negative tracks, principally due to the net proton 
yield), and a peak at higher (dE/dx) with negative am- 
plitude (net pions). In order to extract the particle yields, 
Gaussian functions were fitted to the two extrema in each 
difference distribution. A correction to the net proton 
yield for the difference between K + and K~ yields, based 
upon NA49 data |jl9), was 15% at y cm < and was neg- 
ligible for y cm > 1 . The negative hadron yield was deter- 
mined from the yield of all tracks from negatively charged 
particles excluding electrons, identified via (dE/dx). 

Corrections for detector acceptance and track recon- 
struction efficiency were calculated by embedding and 



reconstructing simulated tracks in real events |18[. The 
tracking efficiency was greater than 95% over most of the 
acceptance, falling below 80% in regions of high track 
density and near the edges of the acceptance. Yield as a 
function of rapidity was determined only for bins having 
acceptance to px=0, with the exception of the two lowest 
rapidity bins for h~ where extrapolation was based upon 
the adjacent bin having full coverage. At high px, the 
acceptance was limited to the region in which detector 
effects were well understood in this analysis, extending at 
least to pt=2-0 GeV/c for all rapidity bins. Instrumental 
background due to secondary interactions with detector 
material was estimated to be 5% of the total measured 
yield, using a GEANT-based Monte Carlo simulation and 
the VENUS event generator jlo| , which reproduces trans- 
verse and forward energy distributions in central Pb+Pb 
collisions ||. Excluding the decay corrections discussed 
below, the systematic errors of all rapidity distributions 
were less than 10%. 

The measured yields contain contributions from the 
products of weak decays that were incorrectly recon- 
structed as primary vertex tracks. The background cor- 
rection to the h~ yield due to the decay of Kg was es- 
timated using NA49 data to be 5% at midrapidity, 
decreasing strongly at higher rapidity. Background to the 
h~ yield due to A decay is less than 1%. The background 
correction to the p—p yield is due to the decays of A, A, 
E + , and £~, and was assessed using a GEANT-based 
simulation of the decay of A and A and reconstruction of 
their charged decay products. (A and A should be under- 
stood in this paper to include the contributions of S° or 
S° as well as feeddown from weak decays, which are not 
distinguishable from primary A or A in this analysis.) To 
investigate the dependence of this correction on the phase 
space distribution of the decaying particles, three differ- 
ent shapes of rapidity distribution for A — A were used, 
based upon (I) predictions of the RQMD model || and 
(II) the VENUS model (h}, and (III) a preliminary NA49 
measurement [^o| . In each case the S + — E~ yield was as- 
sumed to be 30% of the A— A yield |2l]] and to have the 
same rapidity distribution. Strangeness is not conserved 
if both the hyperon yields in I or II and the measured 
kaon distribution are used. The hyperon yields in these 
cases were therefore scaled by the ratio of total charged 
kaon yields in the data and the model, thereby imposing 
strangeness conservation. In order to compare correc- 
tions due to A— A distributions differing only in shape 
but not total yield, the distribution in III was scaled so 
that the ratio of A+S° to charged kaon total yield agreed 
with that from the models. The p—p distribution cor- 
rected by other A — A distributions can be derived from 
the figure. 

Fig. [l], upper panel, shows the event normalized net 
proton yield as a function of rapidity for central Pb+Pb 
collisions, incorporating the three A— A corrections and 
corresponding A — A rapidity distributions. Also shown 
is the proton rapidity distribution for p+p collisions at 
400 GeV jlq] , which is qualitatively different. The com- 
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parison highlights the importance of multiple collisions 
to baryon stopping in nuclear collisions. 

To determine the net baryon (B — B) rapidity distri- 
bution, the contribution of the remaining net baryons, in 
addition to p— p and A — A, must be estimated. Model 
calculations ||[l0| indicate that the rapidity distribution 
of net neutrons follows that of net protons over most of 
phase space, with a 7% larger yield. The contributions 
of £ + — E~ and £~ — S + were accounted for by assum- 
ing the same rapidity distributions as A— A, and scaling 
the A— A distribution by an empirical factor 1.6±0.1 de- 
rived from hadronic reactions pl[ . Multistrange baryons 
contribute less than 2% to B — B ^3|, and were not ac- 
counted for. The net baryon yield B — B was then calcu- 
lated as 

B-B = (2.07±0.05) • (p-p) + (1.6±0.1) • (A-A). (1) 
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FIG. 1. Upper: normalized rapidity distributions of p—p 
for Pb+Pb collisions incorporating correction I (open circles 
are data reflected about y cm =0, errors not shown). Lines 
show variation in data using corrections II (dotted) and III 
(dashed). Also shown are corresponding A — A rapidity distri- 
butions (I solid, II dotted, III dashed), and the scaled proton 
distribution for p+p collisions. Lower: normalized rapidity 
distributions of B — B from Eq. [j] for Pb+Pb incorporating 
correction I, and scaled B — B for S+S. Lines correspond to 
corrections II (dotted) and III (dashed). 

Fig. 0, lower panel, shows the B — B rapidity distri- 
bution for Pb+Pb. The data points and errors are for 
A— A correction I. The variation in B — B due to the use 
of different A — A corrections is smaller than that due to 
other sources of systematic error. The B — B yield is less 
sensitive than the p—p yield to the assumed A— A distri- 
bution, because the A— A distribution is removed from 
the measured p—p distribution by the decay correction 



but added again in Eq. Qto calculate B — B, with the two 
contributions approximately cancelling. Conclusions on 
baryon stopping are therefore based upon the net baryon 
rather than net proton rapidity distribution. Also shown 
in Fig. 0, lower panel, is the B — B rapidity distribution 
for the 3% most central S+S collisions at 200 GeV per 
nucleon | fi2[ , using a coefficient of 2.0 for p—p in Eq. |l| 
and with integral normalized to the number of nucleon 
participants in Pb+Pb. Within \y cm \ < 2.5, there are 
352 ± 12 participant baryons for Pb+Pb and 52 + 3 for 
S+S central collisions p"2| ]. 

The B — B rapidity distribution is narrower for Pb+Pb 
than for S+S collisions, indicating increased baryon stop- 
ping for Pb+Pb collisions. (In the CM frame, ybeam = 2.9 
for Pb+Pb and 3.0 for S+S.) The median rapidity shift of 
the leading nucleon with respect to the beam in high en- 
ergy p+Pb collisions is 2-2.5 units (22|. Due to the sym- 
metry of the nuclear reactions studied in this work, dis- 
tributions of target and projectile nucleon rapidity shifts 
in full phase space cannot be determined separately, lim- 
iting the magnitude of a calculated rapidity shift even 
in the case of large baryon stopping. The mean rapid- 
ity shift relative to ybeam for participant baryons within 
< y cm < 2.5 is (Ay) = -1.76 ± 0.05 for Pb+Pb and 
(Ay) = -1.63 ± 0.16 for S+S collisions. 
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FIG. 2. Normalized rapidity distributions of h~ for cen- 
tral Pb+Pb collisions, and scaled central S+S and isoscalar 
inelastic N+N collisions. 

Fig. ^| shows the event normalized rapidity distribu- 
tion (assuming pion mass) of h~ in central Pb+Pb and 
S+S ]l^] collisions, and isoscalar inelastic N+N [Q colli- 
sions. An independent analysis of Pb+Pb collisions using 
the NA49 Vertex TPCs found a distribution consistent 
with these data Q. The S+S and N+N distributions 
were scaled relative to that for Pb+Pb by the ratio of 
the number of participant nucleons and a factor 0.96 to 
account for the energy dependence of average multiplic- 
ities measured in hadronic collisions p5| . No correction 
for the net isospin difference between Pb+Pb and the 
other systems was applied. Under the assumption that 
the net isospin in the final state of a Pb+Pb collision is 
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carried entirely by mesons, the additional scaling factor 
needed for comparison of the S+S and N+N distributions 
to that of Pb+Pb is 1.12. 

The enhancement of meson yield at midrapidity for nu- 
clear collisions relative to N+N collisions has been noted 
previously [pLlJ| . Taking into account the net isospin dif- 
ference, the agreement between the Pb+Pb and scaled 
S+S distributions is striking. This scaling of the yield 
of produced particles with number of participants, first 
observed in p+ A collisions |2y] , is also observed in the col- 
lision of very heavy ions. Extrapolation of the h~ yield 
to full phase space was performed by taking into account 
the asymmetry in the h~ rapidity distribution (assuming 
pion mass) due to the contribution of K~ jL9| , resulting 
in an estimated total hr yield of 695 ± 30 particles. The 
number of hr per participant nucleon pair, not adjusted 
for isospin, is 4.0 + 0.2 for Pb+Pb collisions, compared to 
3.6 ± 0.2 for S+S collisions and 3.22 ± 0.06 for isoscalar 
N+N collisions lO . 
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FIG. 3. Spectra for central Pb+Pb collisions. Upper: 
transverse mass spectra near t/ cm =0. Solid lines indicate fits 
discussed in text; dotted lines are extrapolations of fit func- 
tion. Lower: rapidity dependence of (px) for hr and p— p. 

Fig. |^, upper panel, shows transverse mass spectra 
(rriT = \J m? + px 2 ) near midrapidity for hr and p— p 
for central Pb+Pb collisions. NA49 has previously re- 
ported the fit of an expanding hadronic source model |2^] 
to midrapidity hr and deuteron mx spectra and hr cor- 
relation functions J(|, giving a freezeout temperature of 
T = (120 ± 12) MeV and transverse expansion velocity 
0± = (0.55 ±0.12). Fig. |, upper panel, shows the fit of 
this model with fixed /3j_ =0.55 to the hr and p— p spec- 
tra reported here. The resulting freezeout temperatures 
are T= (126 ±2) MeV for hr and T= (118 ±5) MeV for 
p-p. 

Fig. ^J, lower panel, shows the rapidity dependence of 
(px) calculated within < px < 2.5 GeV. For p—p for 



central Pb+Pb collisions, (p T ) = (825 + 37) MeV at y cm = 
and (600 ± 17) MeV at y cm = 2.4, compared to the 
rapidity-averaged value (px) = (622+26) MeV for central 
S+S collisions jQ]]. For ]iT for central Pb+Pb collisions, 
(p T ) = (385±18) MeV at y cm = 0, compatible with (377+ 
4) MeV for central S+S collisions p2[ . (In considering 
(pt) for h~ at high rapidity, note that the kinematic limit 
for pions for N+N collisions at 158 GeV falls below px — 1 
GeV near beam rapidity). 

A similar characterization of px distributions re- 
sults from fitting (1/mx) ■ dn/dmx with a function 
A ■ exp(-m T /T). A fit within < to t -to < 0.8 GeV 
to the p—p data near midrapidity for central Pb+Pb col- 
lisions yielded T = (308 ± 15) (NA44 reported (289 ± 7) 
MeV for protons |l3[). NA35 found rapidity -averaged 
T= (235 ± 9) MeV for central S+S collisions Q. The 
increase in T and (px) with particle mass is consistent 
with a larger transverse radial flow velocity at midrapid- 
ity in Pb+Pb than S+S collisions §jl|. 

In summary, baryon stopping for central collisions at 
ultrarelativistic energies increases with system size. The 
integrated yield and rapidity density of negative hadrons 
exhibit scaling with the number of participant nucleons 
for nuclear collisions, and a small enhancement with re- 
spect to N+N collisions. At midrapidity, a large increase 
is observed in {px) for the stopped baryons in Pb+Pb 
relative to S+S collisions, with no significant increase in 
(px) for negative hadrons. The increase in (px) with 
particle mass is consistent with an increase in transverse 
radial flow velocity for heavier colliding systems. 
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